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In the first part of the papei we calculate and discuss the distribution function of the 
electrons of a slightly ionized plasma under the influence of an external magnetic held 
and of an e.m. wave of the type E=E cosooL 

It is shown that, taking into account both elastic and inelastic collisions between elec- 
trons and molecules, it is possible to calculate explicitly the mean electronic energy due only 
to the absorption of the extraordinary wave of the electric field and that it is maximum at 

ed~t 

the gyromagnetic resonance (co = co#= — — ). The study of the effects of inelastic collisions 

shows that in order to have an equal electronic temperature, the intensity of the electric 
field must be, in the case of inelastic collisions, about five times greater than that for elastic 
collisions only. 

In the second part of the paper we calculate, for the wave E = E cosco£, the components 
of the complex dielectric permittivity tensor, that are given by expressions dependent in a 
rather complicated way on E . It is shown that if Eo is sufficiently small they can be simpli- 
fied and, using a suitable perturbation technique, it is possible to calculate a nonlinear dis- 
persion relation. This relation contains explicitly E , and, for E — >(), becomes the Appleton- 
Hartree formula. Furthermore, it shows that in the nonlinear case, too, the electric field 
is split up into two components whose propagation is never independent. 

In the third part of the paper we calculate and study the electronic distribution function 
and the complex dielectric permittivity tensor for a wave of the type E = E () [1 + 77COs(a£-f-(3)] 
COSw£ with co^o>ff»a. 

1. Introduction 

The study of nonlinear propagation of electromagnetic waves in plasma has recently 

aroused the interest of many scientists. In general, the methods used in describing the process 

of propagation of electromagnetic waves in plasma are based on a procedure of linearization 

of the equations that give the mathematical representation of the physical problem. A very 

important physical characteristic of the aforesaid nonlinear effects rests on the fact that they 

can he produced by relatively small electric fields. This can be understood quite easily if we 

consider a slightly ionized plasma in which the electron-electron and electron-ion collision 

frequency is much less than the electron-molecule collision frequency. In fact, the propagation 

of an electromagnetic wave in a plasma causes a relatively high increase of the kinetic energy 

of the electrons. This is mainly due to two facts: the first is the quite large value of the mean 

free path X of the electrons in the plasma (so that they can acquire considerable energy from 

the wave between two collisions), the second is that due to the smallness of the ratio 

2m 
5=- TT c^3A • 10~ 5 between the mass of the electron and the mean mass of the molecules, the mean 
M 

energy transferred in a collision from the electrons to the neutral component of the plasma 
is almost negligible. As a final result we see that the energy distribution of the electrons 
is altered (with an increase of their temperature and mean energy). This variation depends 
on the quantities E, H, co that characterize the wave in the plasma. So that the parameters 
that characterize the plasma from the electromagnetic viewpoint (as the dielectric permittivity 
e ik , the susceptibility xuc, the conductivity a ik and so on) will depend on these parameters 
(E and co mainly) and on the parameters that characterize the plasma from the kinetic view- 
point (as the mean free path X, the collision frequency v and so on). In this way in the rela- 
tions that give the electric polarizability P or the conduction current density j c that is 

Pt=x ik E k 

j c . = a ik E /c 
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the electrical parameters e ik , xuc, and o ik will depend on E also and the quantities P t and j e . 
will no longer be proportional to E. Therefore, the electrodynamic processes in the plasma 
and the same propagation of electromagnetic waves will become nonlinear. There are two 
types of approach to the theory of nonlinear effects in plasma: the first one is a generalization 
of the mean free path method and the second one is the statistical method based on Boltzmann 
equation. The fundamental equation of the first method is Langevin's equation 



- m [ E+y e X ^]' M 



Tt +vy 



e ) m, v, v being respectively charge, mass, collision frequency, and velocity of the electron; E 
is the electric field of the radiowave, ${ the external magnetic field in which the plasma is 
immersed. Knowing E, (1) allows us to calculate v, furthermore putting by definition: 

}t=nev 

where n is the electronic density and j t is the total current density, it is possible to calculate 
the conductivity o ik and the dielectric susceptibility xm by means of the equation 

In this way are obtained the expressions of the dielectric permittivity e; A .(e a -=47r[xa + <5a]) and 
of the conductivity a ik that must be substituted in the equations of the wave propagation. It 
is easy to see that these expressions of e ik and of a ik do not depend on E and so the propagation, 
also of a very strong wave, is linear. To overcome this difficulty [Ginzburg and Gurevich, 1960], 
the electronic collision frequency v (that enters into the expressions of e ik and of <j ik ) is considered 
dependent on the electronic temperature T e , e.g., in the following way, 1 



Kr)=F<»^5? 



T 

being T the molecular temperature and z> (0) the collision frequency in the absence of an electric 
field (T e =T). In order to calculate T e a differential equation is established equating dTJdt 
to the net power gained by the electron, in the following way: 

Given the rather complicated structure of this equation, the usual method to solve it consists 
in disregarding, at first, the dependence of v on T € . In the final formulas obtained in this way 
is introduced v, given by (2), and so it is possible to calculate T e as a function of the various 
parameters of the plasma among which is the electric field. In this way e ik and a ik depend on 
E and the equations for the propagation become nonlinear. It is clear that this type of approach 
to the theory of nonlinear propagation presents defects and contradictions that are typical of 
the mean free path method. This method, worked out particularly by Townsend and Huxley, 
is based on the calculation of the distance traveled by the electron under the influence of the 
accelerating field during the interval of time between two successive collisions. The mean 
distance traveled by an ensemble of electrons per unit time gives the diffusion velocity of the 
electrons through the gas. The intensity of the field is always supposed to be weak: the calcu- 
lation is done by supposing that all the electrons have the same velocity, then the mean value 
of the results is calculated by means of an electronic distribution function, previously unknown. 
A maxwellian distribution function is employed to obtain results which do not contain the above 
mentioned mean values. This method has the drawback of being valid only for weak fields, 



i This dependence is valid for electron-molecule collisions [Ginzburg and Gurevich, I960]. 
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which are not able to deviate the electronic distribution function from its stationary form 
(maxwellian), without the possibility of exactly defining the general validity and the order of 
magnitude of the employed approximations. Furthermore, according to Huxley [1937, 1940, 
1951, and Huxley and Crompton, 1962,] the method of the mean free paths can be employed 
whenever we suppose that among all the various electron-molecule types of collision the binary 
ones are by far the most important and that the motion of the particles can be divided into 
short periods, in which the collisions occur, separated by comparatively long intervals during 
which the interactions with the other particles can be neglected in comparison to the action of 
the field. It is evident that if the behavior of a physical system can be outlined in the above 
mentioned way, it can be analyzed by means of Boltzmann's integrodifferential equation. 
From a more general viewpoint, we do not see how it would be possible to base the theory of 
nonlinear propagation on the equation (1) of dynamics, when we are faced with an electron gas, 
with a statistical velocity distribution. 

1.1. Fundamental Equations of the Statistical Method 

We are interested in the lower flayer of the ionosphere, i.e., the zone between 80 and 95 
km height. The molecular temperature of this layer varies slowly from 205 °K (at 80 km) 
to 217 °K (at 95 km). Let us suppose that the temperature gradient is zero and let us take a 
mean molecular temperature equal to 210 °K. The molecular concentration is of about 10 14 
niol/cm 3 in this layer. The electron density varies from 5.10 2 el/cm 3 (at 80 km) to 5. 10 4 el/cm 3 
(at 95 km) and therefore it is much inferior to the molecular one. The statistical state of the 
ionospheric plasma is described by means of two distribution functions: one /(r,v,/) for the 
electrons, the other F(R,Y,t) for the molecules; these are the solutions of the following system: 



*1 
dt 



+ WrJ+yVv.f=Cn+C V2 



■ ~ + V • VR^+r • VvF=C n +Cn 

in which C n represents the collision term for electron-electron interaction, C 22 the collision term 
for molecule-molecule interactions and last C i2 and C 2l , the collision terms for electron-molecule 
interactions and where we have indicated with y and r, the external accelerations acting on the 
electronic and molecular gases, respectively. Taking into account the characteristics of the 
ionospheric plasma, the second of the written equations can be simplified in the following way: 

C 22 = 

which, after integration, gives the following (maxwellian) distribution function for the molecular 
velocities. 



r ~ iy \2wkTj exp \ 2kT J 



where N is the molecular density. Moreover, as N is much greater than n (electronic density), 
the first equation of the system can be simplified, neglecting C n with respect o C 12 ; in this way 
we obtain: 



\+y • Vr./'+^ [e+ J X ^]-Vv/=tf,2 (4) 

C !2 is given by the following integral operator that we will indicate with J{f] 

c»-J{ J }= JJ[/(v')/ / (V')-/(v)/''(V)i r /«r( ( ? , g )dVda (5) 

where g— |v— V| is the modulus of the relative velocity, <r(#,g) is the differential cross section, 
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# the scattering angle, v' and V are the velocities of the electron and of the molecule before the 
collision (after the collision they become v and V respectively) . Taking into account the char- 
acteristics of the ionospheric plasma, it is possible to simplify (4), [Davydov, 1936 and 1937]. 
We have already said that <5, the mean energy lost by an electron in a collision, is much less than 
1. In the case of elastic collisions, 

•-£=*.. 10- 

in the case of inelastic collisions the energy lost by collision is noticeable, but because these 
collisions are much less frequent than the first ones, the mean value of 8 is always much less 
than 1. For this reason, the rms velocity of an electron is much greater, also in strong fields, 
than the mean directed velocity, so that the symmetric portion in v, of the electronic distribution 
function is much greater than its antisymmetric portion. Moreover, as 8 is small, the collision 
term (5) can be expressed in differential form; in this way (4) becomes the Boltzmann-Fokker- 
Planck equation. For these reasons let us develop the electronic distribution function in 
spherical harmonics in velocity space, putting: 

/(r, v, 0=/o(r, v } 0+a, -li(r, v, t)+x(r, v, t) (6) 

v 

where a « == o' If we put this equation into (4) and integrate all over dtt (differential of the solid 

angle in velocity space) we obtain: 

Putting again (6), multiplied by a y into (4) and integrating again all over the solid angle dtt 
we have: 

^T+rVr /o+~ ^° E+^ W. Xf,+rfi=0 (8) 

ot m ov mc 

where we have put: [Ginzburg and Gurevich, 1960; Davydov, 1937; Chapman and Cowling 
I960]: 



W.!=^K>+^>] 



J{fl}=-rfl. 

In this way, starting from (4), with the approximation (6), we have obtained (7) and (8). To 
establish these equations we employed the properties of monodromy and orthogonality of 
/o and fi expressed by (6) and we supposed that the function % and the functions that can be 
obtained by it with the application of the various operators of (4) are very small in respect to 
the function f and of the corresponding ones, obtained by applying to it the various operators 
of (4) . It is easy to see that if the electronic densit}^ varies smoothly and if the variation of 
the electronic current along a mean free path X is small with respect to the product of the 
electronic density and the velocity, that is if, 



«n(v) 



then 

X~V5|fi|^5/ 

so that the function % is small not only with respect to/ (as we supposed to derive (7) and (8)) 
but also with respect to |fi|. It is clear that with the system composed of (7) and (8) it is 
not possible to determine corrections to/ of the order of 5/ . 
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2. Electronic Distribution Function for a Monochromatic Wave 2 
Let us suppose that the plasma is homogeneous, in this case (7) and (8) become, 



dt '3mv 2 dv 



2v 2 dv \ L mv dflj* J 



^-^E+-^ffxW=0. (10) 

ot m dv mc 

In order to solve this system of four equations we must define two relaxation times. The first 
is the relaxation time t e of the electric field that is defined as the time necessary to change 
substantially the field (if E=E cos ot then t^/^co -1 ). The second is the relaxation time r r of 
the electronic energy; r T is of the order of {bv)~ l as can be seen by integrating (3) with E =0. 
Following Gurevich [1957], we can distinguish two cases: the first in which the electric field 
varies slowly (r^>>r r ), the second in which it varies quickly (r#<<T r ). Now we must 
consider the case in winch the field varies quickly because we study the phenomena related 
with the propagation of a monochromatic wave: 

E=E cosco* (11) 

which pulsation is equal, or very near, to the gyropulsation of the medium. In fact the value 
of the magnetic field of the earth, that enters into (10) varies, in the zone of the ionosphere 
in which we are interested, between 0.35 Gauss and 0.45 Gauss so that, being w^co//, t e c^1A-10~ 7 
sec, because oo H =e | (K \/mc; moreover T r c^(8v)~ l is about 4 -10~ 3 sec, being 5^^3.4 -10~ 5 and vo^S -10 6 

coll/sec. From these calculations we see that in our case ( -- 1<<1. Therefore let us develop 
J and fi into a power series of the parameter (— ) putting, 

jo— /0O+./0I+/02+ • • • 
fl = fl0 + fll+fl2+'... 

we find that in the zero order approximation we can neglect the variation of the distribution 
function due to collisions (of the order of//r r ) in respect to the first term of (9)( because of the 



order of —^JIte ) so -^=0. Therefore 



/oo=JooW 



that is, in the zero order approximation, the symmetric part of the distribution function is 
independent of time. If we put into (9) 

f - u d ^ 00 
tl0 -~ u ^zT 

we see that u satisfies Langevin's equation (1). This equation can be solved for E given by 
(11) and gives the following expression for the persistent part of fi 

fio=A(«) ^° cos »*+B(«) ^° sin «*. (12) 



2 Some results of this section are already contained in the paper "Sull 'eccitazionc dell' "airglow" per mezzo di radioonde", by P. Caldirola, 
Published in Nuovo Cimento Suppl. XIX Series X No. 2 (1961). 
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The two vectors A(co) and B(co) are given by: 

1 



A ( \ e Eo 



.{- 



io H (w 2 H —co 2 -{- v 2 ) 



^XE 



B(«)=- 



I / 2 p 2 1 2\ E l ^ 2 Q 2 I 2\ "J COs2 j »l \ 

+"(«H+« 2 +^) £- + „("„- 3o, 2 +,< 2 ) (0 2 +y 2 i^j y 



e-go 1 

m [(w+w H y+ v 2 ] [(co-co ff ) 2 + 



" 2 ]\ 



-2coco#j> 



W\Eo 



(2 2 2\ E 0| / 2 2 1 Q 2\ ^COS 2 $ <K 



E 



" 2 +" 2 m\ 



}\ 



(13) 



where $ is the angle between Eq and ${ . Given that J (v, t)c^Joo(v)+J i(v, t) and/ o(i')>>/oi (v, 
we have: 

a>jUM)_ ^ & 



d£ 



from which 



i r l 1+ ^ 55; /o ~3^ ^ E -H = ^ * 

\2z> 2 dv J 



a/ 01 (M)=H 2^^- w 



which is not limited for t— > + °°. Stipulating then the existence of the electronic distribution 
function we must impose the condition 3f = 0, that is, 



(,JcTb\ 4 2e 



E.f IO =0. 



This is the differential equation for/ o- Seemingly it is not self-consistent because the time is 
explicitly contained in the expression E-f 10 , but noting that / 00 must contain only the rms 
value of the electric field, we obtain, 



2e 



3md 



(E'fi )T E =V 



e 2 El , . 4/oo 
3^5 ^ W 



having put 
so that, 



t N cos 2 <i> . sin 2 <£ 

<pw = 21 2 H — ^~ 



\(o,+o> H y+v 2+ (o,-o> H y+v>f 



J m =Cexp * 



-r 



mvdv 



e 2 El 



(14) 



(15) 



where the constant C must be determined from the normalization condition, 

/» 00 

I 4irt%o(0)dtf=l. 

We see that the method suggested by Gurevich is particularly clear and precise, because the 
distribution function in the 71th approximation is obtained by imposing the condition of bound- 
edness for t— > + oo of the subsequent approximation. Keeping in mind that the magnetic 
field 3K produces a plasma anisotropy breaking up the plane polarized electromagnetic wave 
into two elliptically polarized waves, it is important to get the electronic distribution function 
for an elliptically polarized electric field E. We express E in terms of its components along the 
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three principal polarization axes, 

E=E, |0 cos a>t + Ei +e iat + EZoe- i<at 

where E|| is the plane-polarized field parallel to £H and E~ and E+ are two circularly polarized 
fields in a plane perpendicular to ${ and rotating, respectively, in the same sense or in the opposite 
sense of the electrons in the magnetic field. For such an electric field the distribution function 
(15) is still valid provided that Elviy) is replaced by, 



Efo«- 



F 2 



2 +v 2 ' (u + a: H ) 2 + v 2 ' (w-co„) 2 + , 2 



(16) 



It is clear that at the gyromagnetic resonance (coc^co#) the energy transferred to the plasma is 
mainly due to the action of the E^q wave. Therefore, let us consider in E$<p(v) only the part due 
to this wave, in this case (15) becomes, 



/ooO) = <°ooexp < - 

Jo kT- 



= C c 



1+ 



tr 



[(«- 



mvdv 



2e 2 E~ ± l 1 

'3m 8 (oo—u H ) 2 +v 2 \~ 2 j 



*h) 2 -\ 



sbwt 2 



where we have put, 



ZmSkTJ _ 



i \ v 



{mv 2 ^ 
-wtf 



(17) 



(IS) 



with X= (7raoA^) -1 being a the "radius" of the molecule, N the molecular density and X the mean 
free path of the electrons that we take as a constant. According to some authors (Phelps, I960; 
Sen and Wyller, 1960] it would be better to choose v proportional to v 2 . We have chosen a 
linear dependence for the greater simplicity of the final formulas obtained, this is tantamount to 
substituting in them a suitable value X e ff, for the mean free path, that takes into account the 
dependence of X on v. Putting, 



15k 2 T 2 



m\ 2 , 



-»*) , ) 



by means of (17) we can calculate the mean energy <€> absorbed by the electrons 
of the plasma under the action of the component Ej_ of the electric field, we get: 



<e>=i*r(7+/0 



. v/ ^itH 



w (hl'bb^) 



(19) 



where W(k, m; z) is the Whittaker function [Whittaker and Watson, 1958] of parameters k, 
m, and argument z. From (19) we have, 






(19.1) 



so that the energy absorption is maximal at the perfect resonance, w = co# (see figs. 1 and 2). 

744-7 58— 65 4 39 




0.3 0.71 
0.5 0.9 



Figure 1. Values of <e> calculated by means of (19) in 
which: 8 = 3.4-10~ 5 , co# = 7. 888-10^ puis I sec corresponding to 
$-£ = 0.42 gauss, \=10 cm corresponding to a pressure 
p = 6.7-10~* mmHg, T = 2W °K. 

It is possible to see clearly that the energy absorption maximal for u=coh. More- 
over for weak electric fields and for increasing — all the curves tend to the same 

OJH 

horizontal asymptote that is equal to the value of the energy of the electrons in 
thermal equilibrium with the molecules. The various curves conespond to the 
following values of the i?_iTo component of the electric field: 

Curve 1 E-l. o=10- 3 volt /cm 

Curve 2 #_iio=7.5-10- 4 volt/cm 

Curve 3 J£jro=5-l(H volt/cm 

Curve 4 #_co=2.5-10-4 volt/cm 

Curve 5 J5j:o=10-4 volt/cm 

Curve 6 #_£-<)= 5 -lO" 5 volt /cm 

Curve 7 Ej:o=10~ 5 volt/cm 



Figure 2. Values of <e> calculated by means of (19) with 
the same values of the parameters as for figure 1 . 

For this figure we can make the same observations as those that we made for 
figure 1. The various curves correspond to different values of — namely: 



Curve 1 


03 = UH 




Curve 2 


CO = 0.7 OJH, 


a) = 1.3 ojh 


Curve 3 


CO = 0.5 OJH, 


jd = 1.5 OJH 


Curve 4 


d) = 2 OJH 




Curve 5 


o; = 3 ojh 




Curve 6 


C0 = 4 OOH 




Curve 7 


co = 5 ojh 
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From the asymptotic representation of the Whit taker functions, it is possible to see that 



lim < e ( w )>=lim <a(a,„)> =-;- kT. 



From which we see that if the frequency of the wave is far from the gyrofrequency, the elec- 
tronic distribution function tends to become maxwellian with a temperature equal to the 
molecular one. At the perfect resonance (17) becomes, 



/cio r) W = a r) exp^ 



r 

rnvdv 



-s. 



TfJ<\ 2tf R±_q\ 



V- (20) 



3m8v 2 

If the electronic temperature is much greater than the molecular one, that is, if 



2 kl< ■ mtf 



(20) becomes: 



ya rt w=w«p{-|<^)V}. (2D 



We have Druyvesteyn's distribution function [Druyvesteyn and Penning, 1940], wherein 
appears the rms value of the electric field as it should be. With the aid of (21) we can calculate 
the fraction P* m (v>v ) of electrons having a velocity exceeding a certain value v B and the 
mean energy <«* <r) > of the electrons. These quantities are given by, 

w ( l r 3.3 / mi? V\ 
P *iru v>v)= _ V 8?8 ' 8 KeEloK) J 

K!) /2 

Measuring <e* (r) > in eV, EZo in volt/cm, X in cm and the pressure p in mmHg, we have: 

< e *(')> = 103.58#I X=6.73^°- (23.1) 

We have found several expressions for the electronic distribution function in a slightly ionized 
gas under the action of an alternating electric field and a constant magnetic field, considering- 
only elastic collision between electrons and gas molecules. As a result of such calculations 

[see fig. 3 and eq (23)] it was found that for — — ~1 (that is, for a field intensity of some milli- 
volt per cm and a pressure of some 10" 3 mmHg) it is possible to have electrons accelerated in 
such a way that their mean energy is several eV so that they can collide melastically with the 
molecules. 

Let us now study the effect of these collisions on the electronic distribution function. 
In a number of researches (Davydov [1937], Druyvesteyn and Penning [1940], Smit [1936]) 
approximate methods have been developed in order to take into account this effect and recently 
a detailed study has been carried out by Kovrizhnykh [I960]. The result achieved by him 
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solves, in principle, the problem raised by us, even if the final formulas are generally too cumber- 
some to be numerically calculated. Due to its laboriousness and to the uncertainty of the 
experimental data required for its application, Kovrizhnykh's method is difficult to be em- 
ployed for practical purposes. Consequently, we thought it advisable to adopt, for the evalua- 
tion of the effect of inelastic collisions on the electronic velocity distribution, a semi-empirical 
procedure particularly useful at least as far as our problem is concerned. Let us consider a c as 
(like air or the ionosphere) consisting essentially of poliatomic molecules. Tn this gas not only 
the levels of the electronic configurations but also the rotation and oscillation levels whose 
energy is rather low (of the order of (10" 4 -10- 2 ) eV for rotation levels and (0 1-0 5) eV for 



1.25, 




Figure 3. Shape ofjhe function P*« (v>v ) given by (22) 
for 8=8.4-10-', E ±Q =W-> volt/cm, p = J0-* mm Hg (that 
is for Ej_ /p— 1 corresvonding, according to (231) to a 
mean electronic energy equal to 6.73 eV), \=66 cm.' 

From this figure it is possible to see that in this case we are compelled to take into 
account also of inelastic collisions. 



in Volt /cm 



Figure 4. Behavior of < € > (for w = co H ) as a function of El,. 

Curve 1: values of <«> calculated by means of (19) in which 5=3 4-10-5 
e urve 2: values of <e> calculated by means of (23.1) in which 5=3 410-5 
L,urve 3: values of <e> calculated by means of (19) in which 5 = 1 6-10-3 
ourve 4: values of <e> calculated by means of (23.2) in which 5=1 6-10-3 
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oscillation levels) can be excited. Therefore inelastic scattering with molecular excitation 
takes place also at low temperature and can absorb a considerable amount of electron energy. 
Keeping in mind the expression of the collision term for elastic collisions, 

we write down in the same way the corresponding term for inelastic collisions, 

where 7?(y) is, corresponding to bv, the electron energy loss due to inelastic collisions. We can 
then take into account this term, by substituting in the formulas for the energy distribution 
of the electrons for elastic collision only, the expression bv with 

B tot (tO=*K«0+«(«0 

which corresponds to the introduction of an "efficient" energy loss for collisions given by, 







(eff) _ ' X (tot) (V) 



v' ew (v) 



The energy loss 5 ieW (v) as a function of the electron velocity has been experimentally 
determined in many papers (for instance [Crompton, Huxley, and Sutton, 1953]). From the 
experimental data one can deduce that, up to energies of about 2 eV, for air and the ionosphere 
b iett) is nearly constant and in both cases it corresponds to about 1.6- 10" 3 (to be compared with 
5 = 3.4 -10 -5 for elastic collisions only). Let us now note that in the formulas for electronic 
energy distribution in the case of elastic collisions only, the parameter 5 appears only through 
the ratio E 2 /b. This leads to the conclusion that, in order to take into account the effect of 
inelastic collisions, it is sufficient to substitute E 2 /8 with E 2 jb {ett) in the 4 final formulas. Thus 
we can say that the effect of the inelastic collisions is to reduce the value of E 2 to a value Z? 2 (red) 
such that E 2 iTed) /b=E 2 lb {ett K We have then 

E (Ted) =E -y^y 
For the air and the ionosphere we shall obtain 



™ „ / 3.4-10- g _l w 



The electronic mean energy is therefore reduced (see fig. 4) from its value (23.1) for 

to 

( € **cr>) 1 (e* ( ^>=20.72#I A=1.345 ^2- (23.2) 

3. Nonlinear Dispersion Relation 3 
Taking (6) into account we have, 

e— 1^ 



J. 



a — loo' 



4t 



E=— en\ v z fiodr. 

4 Jo 



3 A. Airoldi of our Institute has contributed to the work contained in this section. We take the occasion to thank her for the aid given in the 
numerical computations. 
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Putting the z axis along the direction of the external magnetic field a\ we have, 



3 Jo co 2 +v 2 dv 



'-'*»- 3~Jo "2\ 



;,}£•* 



^27/ — 0"l/. 



(co+co^) 2 +j> 2 (co— co^) 2 +^ 2 J dz; 



(oj+co h ) 2 + p 2 (co— co#) 2 
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where co p is the plasma frequency, that is, 



(C0+C0tf) 2 +1 
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(24) 
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The knowledge of the expressions of e a . and <r tt allows us to study the problem of the propa- 
gation of electromagnetic waves in a plasma. In order to solve this problem, it is possible to 
use in the linear approximation two methods that, in general, are equivalent. The first con- 
sists in the integration of the equations of the wave propagation, the second is related to the 
calculation and the study of the dispersion relation. The difference between these two methods 
lies in the fact that: while the first one (the integration of the equations of the wave propagation) 
allows us to obtain quantitatively exact results only in relatively few cases; the second, on the 
contrary, gives a much more general description, also if it is only qualitative, of the propagation. 
More precisely we have that in the case of a uniform plasma the two methods that give equal 
results, in the nonuniform case the knowledge of the dispersion relation allows us, very often, to 
write down approximate solutions of the equations of the propagation. These solutions are, in 
general, the first term of a geometrical optics series. We must stress the fact that if we base the 
analysis of nonlinear propagation of a monochromatic wave launched in a plasma on (24) we 
are compelled to admit that the deviations from linearity are not very strong. In fact we must 
suppose that there is no harmonic generation so that the main effect of nonlinear propagation 
consists in increasing the absorption coefficient and in modifying the velocity of the wave and 
not in a substantial variation of the shape of the monochromatic signal (caused by harmonic 
generation). For these reasons it ought to be possible to write down an approximate solution 
of the wave equation in the following way, 



E=E exp A i-(ct— \ndn V 



(25) 



where now, 



n=n(E ) 



it is clear that for E -^0 this formula must become the Appleton-Hartree formula. 

We are now faced with the problem of calculating the integrals (24). If we consider the 
collision frequency v a linear or a quadratic function of the velocity of the electron it is easy to 
see that, in the case (K ^ 0, it is not possible to give closed expressions for the components of 
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e and a. Furthermore, putting y=const., it is possible to calculate all these integrals that are 
given by expressions that do not contain E () . On the other hand if we evaluate, in this case, 
the mean energy of an electron we easily arrive at the following formula 

so that, also if we have a noticeable heating of the plasma, the propagation is linear. All 
these difficulties can be overcome if we consider only weak fields for which the electronic 
distribution function does not appreciably deviate from its stationary form (maxwellian). 
In this way it is possible to calculate this function by developing /,„, in Taylor series of Ef } , 
stopping the expansion after the first two terms and taking then v= const. If we indicate 
the new distribution function obtained in this way with/(-y, E Q ) we have 



that is: 



j{v, E )=Mv, ^=°)+(|^)^ E ° 

■to *Ks*r) L 1+ i>wd exp \~2kff 



the mean energy of an electron calculated by means of this function is given by, 

<«>=!*r[i+|a] 

where, 

2mbkT 
so that we define as a weak field a field satisfying tlie following inequality, 4 



(26) 
(26.1) 

(27) 
(28) 

(28.1) 



By means of (26.1) it is possible to calculate all the integrals (24) and the following expressions 
are obtained, 
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(24.1) 



4 It must be said that this inequality can be too strong, what is important is that the second term at the right of (27) does not become much 
greater than the first. 
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where: 



A=l+a 



(29) 



for EqC^O, that is A=l, these formulas become those of the linear approximation. Before 
calculating explicitly the dispersion relation it is better to find the expression of the components 
of | and g in the usual reference system for the study of wave propagation. To this aim let us 
put the external magnetic field H in the plane yOz of the new reference system at an angle 9 
with Oz axis. It is easily seen that the matrix that allows us to pass from the old reference 
system to the new one is given by 

1 



cos 6 
-sin 9 



sin 6 
cos 9 



If we indicate with e ih the components of the complex dielectric permittivity tensor (de- 

4:7T \ • V 

fined by e ik =e ik — i — a ik ) in the new reference system and with e' ik the components of the same 
co / 



tensor in the old system, we have 
Introducing for simplicity's sake the notation 






12 = co — %v 



x=A 



we have: 



coft(13 2 



l»H) 



e xv = — e' yx = iaQw H cos 9 
e'xz = — e 2Z =ia£La) H sin 9 
e' yy = l-a(tt 2 -u 2 H sin 2 9) 
e yz = 6zy = — a co % sin 9 cos 9 
e' 2Z = l-a(tt 2 — co 2 tfCOS 2 9). 
Let us consider the two Maxwell equations: 

vxH= ^+T J « 



VXE=- 



c dt 



for the plane electromagnetic wave (25) propagating along the z axis they give, 



5 2 E 
dz 2 ' 



tf = 



E=0 



(30) 



the two Maxwell equations and (30) are equivalent if we suppose that there is not harmonic 
generation of the fundamental frequency co. Furthermore we have, 

dE . - <* T? A 

r- — H-nE=0 
dz c 



dz 2 



+^n 2 E= 
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:— % ~r— E 
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but if the electric field is sufficiently weak the propagation will be almost linear so that it ought 
to be possible to disregard the right member of the last equations; this, on the other hand, 
agrees with the spirit of the perturbation technique that we have adopted. 5 So, 



in this way (30) becomes, 



oz 2 c 2 



Wxx — Tl 2 ]E 0x + e xv E 0jl +e xt E 02 = 



e'y X E Ox +[e vu —n 2 ]E 0y +e yz E o =0 



€ f tx E Ox +€ , gy E Oy -\-€ r gg E O!l =0 

where, now, this is not a linear homogeneous system in the three unknowns E 0x , E Qy , E 0z as it 
happens in the linear case, but it is a more complicated system owing to the presence of El in 
the expressions of e ik . However, in this case, too, it is easily shown that there exists a nonzero 
solution for E 0x , E 0in E 0z provided that the following determinant is equal to zero: 



C-n< 



e™— n 



*zz 



= 



from which we have (compare with (6.1) [Budden, 1961]), 

AX 



n 2 =l- 



1-iZ- 
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(31) 
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or (compare with (2.3) [Gullet, 1963]) 
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(31.1) 



For A=l (linear case) (31) becomes Appleton-Hartree's formula. In this case, as it is known 
[Sen and Wyller, 1960], the results obtained by the statistical theory and by the mean free path 
theory must coincide. The two signs that appear in (31) show that in the nonlinear case too, 
the field is split up into two components whose propagation is never independent, given that 
under the square root there appears the total value of the electric field. We will call these two 
waves ordinary generalized (O.G.) wave and extraordinary generalized (E.G.) wave, identifying 
them with the corresponding ones of Appleton (O.A. and E.A., respectively) in the limit E b -^0 
for every propagation condition. Let us put 



W 



e 2 El 



2mbkTu 2 



5 It is clear that we can disregard the right hand side of the last equation if -— 1<<1; this is the usual condition for the validity of W.K.B. 
approximation. 
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and A becomes 

4-1-luT^^ sin 2 $/' 1 1 \"1 . 9 v 

It is interesting to notice that the electronic density n, that is X, does not appear in (32) so that 
all the results that are obtained for the Appleton-Hartree formula when is varied only n are 
valid also now if we substitute X with AX. The zeros of n 2 are located at the points, 

AX+ iZ= 1 AX+ iZ= 1 ± Y 

and one of the values of n 2 is infinite for, 

AX=(l-iZ) {1 _ iZ)2 _y 2 cog2 e ' 

In this way we see that for 4=1 we obtain the well known results of the linear theory [Budden, 
1961]. Furthermore the two values of n 2 are infinite for, 

Y=±(—l±iZ) if sin <£^0 

Z=±i if cos $5*0 

these poles are also branch points for n if 0=~. It is easily seen from (15) and (26.1) that in 

these points we have J 00 (v) = 1 , J(v, E n ) = oo f for these reasons all the integrals (24) diverge. 
Moreover A= °° and in this case the perturbation technique that we have used is no more valid. 
Given that the dispersion relation (31) depends on co in a much more complicated way than 
Appleton-Hartree's formula, we have studied it using an electronic computer. By means of 
(31) we have calculated the refraction index \x and the absorption index X using the relation, 

n 2 =(n-ix) 2 =M-iN 
from which, 

M=[HVM 2 +A^ 2 +M}] 1/2 



x=[^M 2 +N 2 -M}] 



1/2 



In the numerical calculations we have taken as a constant: 

(1) The molecular temperature T, T=200 °K 

(2) the external magnetic field \&{ \, \cH\ =0.42 Gauss for which co H =7 . 388 • 10 6 puls/sec 

(3) the electronic density n, n=10 i el/cm 3 

(4) the mean energy lost by an electron in a collision, 5= 1.6* 10 -3 

(5) the direction of the propagation, we have put = (longitudinal propagation) in this case 

we have, for symmetry, <£= - 

We have varied: 

(1) The amplitude E of the electric field 

(2) its pulsation around the chosen value of oo H 

(3) the electronic collision frequency v. 

The figures show the results that we have obtained. 

From the figures 5, 6, and 7 it is possible to see that the absorption index is in the non- 
linear case greater than that calculated by means of Appleton-Hartree's formula, moreover, 
it increases with increasing electric field. For a field equal to, or less than, 10~ 4 volt/cm we 
do not have appreciable deviations from the linearity. For low values of the collision fre- 
quency (fig. 5, */=2.5-10 6 coll/sec) the absorption index of the O.G. wave increases very much 

for co~co H , with increasing v this increase disappears (figs. 6 and 7) and for equal E and — 

the ratio between the values of % for the O.G. wave (and E.G. wave) and the values of x 
for the O.A. wave (and E.A. wave) tends to 1 with increasing v. 
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Figure 5. Values of x for v = 2.5-10* coll I sec. 

Curves 1 (O.G. and E.G.) # =10-3 volt/cm ai max =2.49 
Curves 2 (O.G. and E.G.) # =7.5-10-» volt/cm a 2niax =1.40 
Curves 3 (O.G. and E.G.) # =5-10-* volt/cm a 3ma x=0.62 
Curves 4 (O.G. and E.G.) # =2.5-10-4 volt/cm ai max =0A5. 



0.5 0.6 0.7 0.8 0.9 



Figure 6. Values of x for v=5-10 6 coll sec. 

Symbols same as in figure 5 

aimax = 0.67 O2max = 0.38 O3max=0.17. 



As it is seen from figures 8, 9, and 10 the behavior of the refraction index /z is similar to that 
illustrated for x when E , w, v are varied. 

The fact that the propagation tends to become linear with increasing v is easily under- 
standable. In fact if E and co remain fixed the power transferred from the wave to the elec- 
trons is unaltered, while with increasing v also the power dissipated by the electrons increases 
due to the collisions, in this way the electronic temperature tends to become equal to the 
molecular one. 
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Figure 7. Values of x for v— 7.5-1CP coll/sec. 
Symbols same as in figure 5 

aimax = 0.32 a 2 max = 0.18 O3max = 0.08. 
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Figure 8. Values of n for v=2.5-10 fS colli sec. 
Symbols same as in figure 5. 

a, ma x = 2.49 C 2 max = 1.40 O 3 max = 0.62 « 4max = 0.15 



From figures 8 and 9 it can be seen that with increasing nonlinearity also the wave number 
of the O.G. wave increases. This fact generalizes a preceding result [Epstein, 1962] according 
to which with increasing nonlinearity, the wave number of an electromagnetic wave passing- 
through an isotropic plasma is also increased. 

The value of a, given by (28) varies along the curves traced in the figures, for every curve 
the maximum value of this quantity is indicated that, in general, is achieved by a for co ^co H . 
It is seen that for some propagation cases the condition (28.1) is not fulfilled. For this reason 
it is interesting to see if we obtain different results taking into account other terms in the 
development of J 00 in series of E 2 Q . To this aim we have calculated the third term of the ex- 
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pansion that gives for A the value 



and the fourth term for which 
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Figure 9. Values of p for v=5-10* coll I sec. 

Symbols same as in figure 5. 
a im.x-0.62 a 2 max=0.38 a 3ma *-0.17 



Figure 10. Values of n for v=7.5-lCP coll/sec. 
Symbols same as in figure 5. 

aimax = 0.32 2rnax = 0.18 SmaX-0.08 
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If we consider the method used to calculate (31) we see immediately that we obtain the new 
dispersion relations taking for A the values given by (33) and (34). The figures 11, 12, 13, 
and 14 show the results obtained in this way. For electric fields equal to, or less than, 10~ 4 
volt /cm, the results calculated by means of the values of A given by (29), (33), and (34) are 
practically equal. For stronger electric fields, of the order of 10" 3 volt/cm, or for low collision 
frequencies (fig. 11 and 13) the results are quite different. With increasing v (figs. 12 and 14) 
these differences disappear, given that the nonlinearity of the medium is decreased for the 
reasons explained before. 

The problem that arises when we want to treat strong fields and low values of the collision 
frequency does not consist in choosing a particular form for A but in showing that, in this 
case, we can still use the dispersion relation. In fact, as it can be seen from figure 8 it is possible 
to find some propagation conditions for which there is a strong difference between the phase 
velocities of the O.G. and E.G. waves. This causes, of course, a strong deformation of the 
monochromatic wave launched in the plasma. In this case we have harmonic generation and 
(25) and (30) loose their meaning. From this analysis we infer that it is not possible to treat, 
using a dispersion relation, the problem of nonlinear propagation of a rather strong electro- 
magnetic wave in a plasma with a rather low value of the electronic collision frequency. We 
think that this problem could be solved in a satisfactory way by trying to integrate directly 
the system constituted by the equations of Boltzmann and of Maxwell. In fact, it is clear, 
that the electronic distribution function (15) from which we started in order to obtain (31) 
is valid only if the wave remains monochromatic. 

4. Electronic Distribution Function for an Amplitude Modulated Electric Field 6 

The system of equations that we must integrate is given by (9) and (10), that is, 

e & r, 2E . fi} s dr 3 / 1+ ^s\ /o i (9) 



dt^Smv 2 dv l 1; 2v 2 dv 



where now 



^ + ^^E + ^^ r Xf 1 +,f 1 =0 (10) 

ot m ov mc 

E=E [1 + 77 cos (at+0)] cos at (35) 



with co#~co>>a, 0<?7<1, 0</3<27r. We begin by defining three relaxation times: the first 
is the quick relaxation time of the electric field r E . q . , it is determined by the frequency of the 
carrier and is of the order of or 1 , that is, te.q.^w" 1 ; the second is the slow relaxation time of 
the electric field r Es ., it is determined by the frequency of the modulating signal, and is of the 
order of a -1 , that is, t^.b.—oT 1 ; the third is the relaxation time of the electronic energy r r that, 
as we know, is of the order of (5j>) -1 . We shall suppose that 

— ^7- E . q .<<7><<>E.8. (36) 

If we remember the calculations of section 2 we see that (36) can be satisfied for values of a 
till to some kc/s. From the fact that r E . s .>>r r derives that now the electronic distribution 
function must depend explicitly on time by means of the function cos (atf+j8). In this way 

we have the first term of (9) of the order of / /-~— — • f Q . Let us put, as for the monochromatic 

* ol r E . s . 

wave : 7 

f- u d/ °- 



6 All the results contained in this section will appear, in a work by one of us (O. D. B.), in a more complete and elaborate form. 

7 It is easy to see that this method of integration is tantamount to developing the functions /o and f 1 in a double series of powers of the pa- 
rameters ( E q ' j «1 and ( E q " ] <<1. We shall calculate the first terms of these developments. 
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Figure 11. Values of x for v = 2.5-10 % coll I sec. 

1(0. G. and E.G.) a with A given by (29) # = 10-3 volt/cm 

1(0. G. and E.G.) a* with A given by (33) 7^ = 10-3 volt/cm 

1(0. G. and E.G.) tf with A given by (34) 77 =10-3 volt/cm 

3(0. G. and E.G.) a with A given by (29) E =5-10-i volt/cm 

3(0. G. and E.G.) a* with A given by (33) # =5-l(H volt/cm 




Figure 12. Values of x for i> = 5-10 C) coll/sec. 
Symbols same as in figure 11. 



We have 






5ud/o u d 2 / 
bt dv btbv 



but ?/ (which gives the current) must vary in time at the rate of the frequency co of the carrier, 
so du/d/^coii, in this way, disregarding terms of the order of a/ a, we have, 



Mi 



dud/o 
dt dv 



and (10) becomes the usual Langevin's equation, 

dt m\_ c J 
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Figure 13. Values of n for v=2.5-10* coll/sec. 

Symbols same as in figure 11. 
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Figure 14. Values of n for v = 5-10 (i coll/sec. 
Symbols same as in figure 11. 



We already know (see sec. (2)) the solution of this equation for a monochromatic wave, taking 
into account that (35) can be written as, 

E=EoJt7 cos {(o)+a)t+^}+E Q COS cot + E %r) COS {(co— a)t — (3} 
and remembering (13) we have the persistent part of fi given by: 

f 1 = [A(a ) )+ir 7 { A(co+a)+A(co-a) } COS (at+0) 

+i rj{ B(co+a) -B(co-a) } sin (a*+|8)] ^ COS ut 

+["B(«)+|i|{B(«+a)+B(«--a)}cos(a*+i8) 

+i7 7 {A( w -a)-A(a;+a)}sinM+/3)J^ sinco^ 
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(37) 



Id the calculation of E • fi we disregard terms of the order of — =-^j that is, we put, 

w T E .s. 

E -A(co-a) = E -A(a;) = E () -A(co+a:) 1 

E -B(co-a) = E -B(a;) = E -B(co+a) J 

and we define a function S(t) by: 

H(0 = [l+^cos(^+/3)] 2 
in this way (9) becomes, 



(38) 



(39) 



- 2 ~ \v 2 { Eo-A(o)) cos 2 co* + E .B(co) cos ut sin a>t }E(0 l^l-A^ [~^ 3 6 +~ ^) /ol=0- 



On account of (36) we see that the first term of this equation (of the order of /o/^e.s.) can be 
disregarded with respect to the second term (of the order of / /r E . q .) and to the third term 
(of the order of / /r r ) so that taking the mean, over the period of the carrier, of the resulting 
equation we obtain, 

£![*{>0££*w»)!}*;h 



a first integration inves 



w%+vv* (kT+ e ^? <p(v)m) |?=A(*) 



where A(t) is an unknown function of time, it can be easily seen that, unless it is not identically 
zero,/ diverges strongly for v=0, so that A (0—0 and 



/ o (0, *)=<?(*) exp h 



-/: 



mvdv 



p 2 E 2 

*r+!#; ^,(„)H(o 



3 w 6 



v. 



(40) 



The study of this distribution function is made easy by the fact that it is identical with (15) 
provided that we substitute El with El&{t). In particular we have the mean energy (e) 
depending on time, (e(t)), s and taking as before only the wave E^ it is given by (19) 
where now 

_e 2 ElW „, 



35FT 2 " { } 



furthermore we have, as before, 






=0. 



Taking into account (38) , (37) for f x can be simplified in the following way 



f i ^ [ 1 + 77 cos (at-{-j3)] [A(w) cos cot-{-B(w) sin oit\ 



dv 



furthermore putting, 



-~ en I *>% 



,^=(T-E + 



dP 



(37.1) 



it is possible to evaluate the conductivity and the dielectric permittivity of the plasma. It is 



' See figure 15. 
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easy to see that disregarding terms of the order of 



Oi t e . ( 



the components of these tensors 



co r E . 5 

are given by (24) provided that we substitute E% with E 2 a{t). We can still use approxi- 
mation (26), so that the components of the complex dielectric permittivity tensor can be written 

as 

€i k (co,a, t) = e' ik (a))+A ik (u,a, t) (41) 

being 

a / i\ > ( \ e 2 E 2 <p(v) . 
A ik {a>,a,t) = e ik (a>) 2 m^T ^ 



where e^(co) is the expression that we obtain in the linear theory for the carrier. By (41) we 
see that e ik is varied, around its steady value e^(co), by the quantity A^(co, a, t). Let us now 
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Figure 15. Behavior of <e(t)>/or E_l c = 

5-10~* volt I cm. 
Other constants same as in figure 1, for various values of 17. 



0000 



OQOOOQOOOOOQQ a ^+B 
v- t- i- CM CM (M rn rn r^, 



56 



hi unci i in the plasma, besides the wave (35), also another wave of the type, 

E r =E 0y cos (yt+A). (42) 

Indicating with <p(co, v) the function given by (14) and with <p(y, v) the same function with co 
replaced by 7 and where <£ is substituted with the angle between E 07 and SK, it is easily seen that 
the electronic distribution function f (v, t) is either: 



r 

f Q (v,t) = C(t) exp h 



i 



mvdv 



kT+^[El<pM3(t) + 2El yV (y, v ) cos 2 (yt + A)] 



dm8 



(43.1) 



or 



/„(», t) = C(t) exp 



,,„ mvdv "J 

q ' 1 "J" kT +^ E fo(»> ")m)+EIMy, v)\ \ (4:5 - 2) 

Equation (43.1) is valid if T E y>y>r n (43.2) is valid if T E <^<^T r , r E being the relaxation time of the 
electric field (42) which is: T E o^y~ l . Furthermore it must be noticed that in the calculation of 
(43.1 and 43.2) we have neglected terms proportional to the product of the two waves (35) 
and (42) ; it is possible to make this approximation if there is no phase correlation between 
the two fields. Let us suppose now that the electric field (42) is small and that 7 is far from 
the gyromagnetic frequency so that we can neglect in (43) the part due to this wave, in this 
way/ (fl, t) is given by (40). The equation for the propagation of the wave (42) is, 

VXVXE 7+ I, d ^ + f%«=0 (44) 

and given that the wave (35) contributes to j ( . ond , the propagation conditions of the wave (42) 
will depend on those of wave (35). We can still use the approximation (26); an easy calcu- 
lation shows that in this case (44) becomes, 

VXVXE T +I l(T ). d ^+f,( 7 ).^4-^=0 

where e(y) and a(y) are the expressions that we obtain in the linear theory for the wave (42) 
and where £f is given by : 






dE 
dt' 



It is clear that the wave (42) does not remain monochromatic and that these considerations 
constitute the starting point for a microscopic theory of cross modulation. 



This paper presents some results of a research program in progress at the Istituto di Scienze 
Fisiche delPUniversita di Milano with the contribution of C. N. R. 
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